Non-perturbative behavior of the quantum phase transition to a nematic Fermi fluid 
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We discuss shape (Pomeranchuk) instabilities of the Fermi surface of a two-dimensional Fermi sys- 
tem using bosonization. We consider in detail the quantum critical behavior of the transition of a 
two dimensional Fermi fluid to a nematic state which breaks spontaneously the rotational invariance 
of the Fermi liquid. We show that higher dimensional bosonization reproduces the quantum critical 
behavior expected from the Hertz-Millis analysis, and verify that this theory has dynamic critical 
exponent z = 3. Going beyond this framework, we study the behavior of the fermion degrees of 
freedom directly, and show that at quantum criticality as well as in the the quantum nematic phase 
(except along a set of measure zero of symmetry-dictated directions) the quasi-particles of the nor- 
mal Fermi liquid are generally wiped out. Instead, they exhibit short ranged spatial correlations that 
decay faster than any power-law, with the law |a;|~'^ exp(— const. |a:|^^^) and we verify explicitely the 
vanishing of the fermion residue utilizing this expression. In contrast, the fermion auto-correlation 
function has the behavior exp(--const. |t|~^'^^). In this regime we also find that, at low fre- 
quency, the single-particle fermion density-of-states behaves as TV* (lu) = N* (0) 4- B uP'l^ log + . • ., 
where N*{0) is larger than the free Fermi value, N{0), and _B is a constant. These results confirm 
the non- Fermi liquid nature of both the quantum critical theory and of the nematic phase. 



I. INTRODUCTION 

The behavior of interacting Fermi systems near con- 
tinuous quantum phase transitions is a central problem 
in the physics of strongly correlated systems. Although 
much work has been done on this subject, there are still 
many open and as yet unresolved questions. At present 
the standard theory of quantum phase transitions^'^'^ re- 
lies primarily on an analysis on the effects of fluctuations 
perturbatively about the results of Hartree-Fock theory. 
This analysis proceeds in almost complete analogy with 
the theory of classical critical phenomena about its up- 
per critical dimension, and its straightforward extension 
to quantum phase transitions. In practice it consists of 
an effective theory for a suitable order parameter field 
while other degrees of freedom, including fermions, are 
often integrated out at the outset. 

In many cases of interest the systems are metallic and 
have gapless fermionic excitations. In the standard ap- 
proach, their net effect is to introduce damping in the col- 
lective modes associated with the order parameter field. 
In practice this results in the introduction of dissipative 
terms in the effective action. While much of this is cer- 
tainly correct, this approach implicitly assumes that the 
fermions are largely unaffected by quantum criticality. 
Why this should be the case is far from obvious. 

The assumptions of the Landau theory of the Fermi 
liquidA^ are self-consistent and well justified within 
the Landau phase which has a sizable basin of stability, 
except in one dimensional systems£i2424iSiiLi2iiiiiiiS4i&. 
However, there is no reason for these assumptions to hold 



outside the Landau phase. However, there is also mount- 
ing evidence that these assumptions may also not hold 
in a number of phases (and not just at quantum critical 
points), including ferromagnetic metals^"* and nematic 
phases of Fermi fluids^^-^". The possibility that quan- 
tum criticality may lead to non- Fermi liquid behavior has 
been a focus of research in recent years, primarily (but 
not only) in connection with the physics of the "normal 
phase" of high temperature superconductorsi2iiS2i2^, 
and with heavy-fermion systemsSi. 

The simplest example where the Landau assumptions 
on the behavior of the quasiparticles are violated is the 
quantum phase transition from a normal (Landau) Fermi 
liquid phase to a nematic Fermi fluid^^. A nematic 
Fermi fluid is a uniform phase of a system of inter- 
acting fermions in which the shape of the Fermi sur- 
face is distorted spontaneously , thus breaking rotational 
invarianceiS^ This state is an example of the fate of a 
Fermi liquid beyond a Pomeranchuk instability^^. In 
this case, the Landau assumptions appear to be violated 
throughout this phase, and not just at the quantum crit- 
ical point liSiSS 

The clearest experimental evidence to date of a ne- 
matic Fermi fluid phase has been found in very clean 
two-dimensional electron gases in magnetic fields in ultra- 
clean samples i2LS£ The striking resistivity anisotropics 
that are observed in these experiments can be explained 
by the onset of nematic order at low temperatures It 
has also been proposed that phases of this type may 
play a central role on the physics of high tempera- 
ture superconductorsi^S*^ This charge-ordered state of 
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a strongly correlated system of fermions is the simplest 
example of an electronic liquid crystal phasa^S. 

The problem of the fate of the fermions at quantum 
criticality, and in the "non-Fermi liquid" phases men- 
tioned above, so far has only been considered within 
perturbative corrections to Hartree-Fock/RPA theory. 
Oganesyan, Kivelson and Fradkin^^ found that the quasi- 
particles are wiped out as well defined quantum states. 
This is due to the large fluctuations of (overdamped) 
quadrupolar collective modes. These authors found, 
within a Hartree-Fock and RPA theory, an overdamped 
collective mode with z = 3 at the critical point. They also 
found that the fermion self-energy acquires, at the quan- 
tum critical point, an imaginary part with a frequency de- 
pendence following the law uP'^'^ . A similar behavior has 
been found in the case of the Stoner transition and in the 
antiferromagnetic phasei^S^. Oganesyan and coworkers 
also found that this behavior holds inside the nematic 
phase, except along a set of measure zero of directions de- 
termined by the symmetry breakingiS However, it seems 
quite likely that such leading order behavio»i2iSli may ac- 
tually signal the complete failure of the Landau theory 
of the Fermi liquid. It is clear that to better understand 
this problem a non-perturbative analysis of the behavior 
of the fermions at the quantum phase transitions (and be- 
yond) is needed. Chubukov'^'^ has given arguments which, 
in the context of the ferromagnetic metallic transition, 
suggest that this behavior may persist beyond the lowest 
order in perturbation theory. 

In this paper we will consider the nematic quan- 
tum phase transition in Fermi fluids using the 
non-perturbative approach of higher dimensional 
bosonization^'''-'^?-^^. We will not discuss the (important) 
lattice effects here. Bosonization is a powerful tool to 
study the non-perturbative behavior of one-dimensional 
gapless Fermi systems, the best understood fermionic 
quantum critical systems^. As it is well known, the 
kinematics of one-dimensional systems is so constrained 
that the bosonic collective modes completely exhaust 
the spectrum of these fermionic systems, allowing 
even for a full reconstruction of the fermionic opera- 
tors entirely in terms of bosons. A striking result if 
one-dimensional system is that the electron acquires a 
non-trivial anomalous dimension and it is no longer the 
quasiparticle of these systems, even for arbitrarily weak 
interactions. For these reasons one-dimensional gapless 
Fermi systems have been termed 'Luttinger liquids'. 
The actual quasiparticles are non-trivial solitons which 
are orthogonal to a bare electron."^^ 

In dimensions higher than one the physics (and the 
kinematics) is quite different than in one dimension. Nev- 
ertheless bosonization methods still yield the physics of 
the Landau theory of the Fermi liquid correctlyiiSiiSii^ 
Superflcially this may seem surprising since in dimensions 
higher than one there are no longer strong kinematic con- 
straints, and consequently the bosonic collective modes 
cannot exhaust the spectrum of an interacting Fermi sys- 
tem. Instead, except for narrow regimes in which the 



collective modes are stable quantum states, they exhibit 
Landau damping, reflecting their decay into particle-hole 
pairs. It is a key check of the validity of higher dimen- 
sional bosonization that it gets the physics of Landau 
damping 

One appealing feature of higher dimensional bosoniza- 
tion is that it is actually a theory of the quantum fluc- 
tuations of the shape of the Fermi surface. It is thus 
a natural approach to study quantum phase transitions 
associated with Pomeranchuk instabilities, and in par- 
ticular the nematic state. ''^ More specifically, we focus 
on the nematic case for spinless fermions and compare 
with the work of Oganesyan and coworkersiS based on 
RPA and Hartree-Fock. We find that the physics of the 
bosonic collective modes is the same in bosonization and 
in RPA, and thus our results agree with those of Ref.jlQlj 
in the Landau phase, in the nematic phase and at the 
quantum critical point. Perhaps, this is not so surprising 
since at long wavelengths RPA becomes asymptotically 
exact and this is the regime in which bosonization is cor- 
rect (for a more thorough discussion, see l39() . In particu- 
lar we derive the effective action near the quantum crit- 
ical point and find that it does have a Hertz-Millis form 
with dynamic critical exponent z = 3, consistent with 
the findings Oganesyan and coworkersi^ and by Nilsson 
and Castro Neto^, but in disagreement with the results 
of K. Yang2&. 

We further use bosonization methods to obtain the 
fermion propagator. This result is well beyond the 
Hartree-Fock/RPA theory and thus it allows us to study 
the fate of the fermions non-perturbatively. We find 
striking violations of the Landau assumptions for Fermi 
liquids. Thus, the equal-time behavior of the fermion 
propagator at the quantum critical point (at zero tem- 
perature) is found to fall off faster than any power, de- 
caying instead with a law jiy exp(— const. |a;|^/'^) as a 
function of distance. The same behavior is found in the 
nematic phase except along symmetry-determined direc- 
tions. We also verify explicitely from this expression the 
vanishing of the fermion residue as expected from this 
kind of behavior. In contrast to the equal-time behavior, 
at quantum criticality the fermion auto-correlation func- 
tion behaves as jiy exp(— const. |f|~^/'^), with a similar 
albeit anisotropic law in the nematic phase as well. We 
also find that the low energy behavior of the one-particle 
density of states N*{io) exhibits an enhancement to a 
zero frequency value 7V*(0) which we find to be larger 
than A^(0), its non- interacting value. At finite but low 
frequency we further find that this one-particle density 
of states behaves as N*{uS) = N*{0) + B w^/s (5 ig 
a constant), i.e. a cusp at a; = 0. 

Thus, our bosonization results confirm that the ne- 
matic phase of a Fermi fiuid is a non-Fermi liquid. How- 
ever, its behavior is more complex than the predictions 
of the Hartree-Fock/RPA theory. Recently Chubukov^'^ 
has examined the behavior of the fermion self-energy in 
perturbation theory at the ferromagnetic quantum criti- 
cal point and found that the frequency-dependence is not 
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changed by higher order corrections. Our results for the 
auto-correlation function are consistent with his results, 
as well as with Refs.^^ and |23|i3|. However we also 
find that the equal-time propagator (the "one-particle 
density matrix" ) has a very different behavior than what 
is predicted from these diagrammatic methods. 

The paper is organized as follows: In section^ we de- 
rive a theory of the nematic QCP via higher dimensional 
bosonization. Here we present a theory of the quantum 
phase transition to the nematic Fermi fluid, section Hi Al 
and show that it reproduces the analog of Hertz-Millis 
theory for this problem. In particular we give a detailed 
analysis of the spectral functions of the collective modes, 
section Hi Bl and derive the effective action valid in the 
vicinity of the quantum phase transition, section FlI CI In 
section IHII we use bosonization methods to calculate the 
fermion propagator. Here we extract the full diagram- 
matic perturbation theory of the fermion Green function 
from bosonization, and use it to calculate the fermion 
self-energy. Here we check that the bosonization for- 
mulas reproduce correctly the non-Fermi liquid behavior 
found within the Hartree-Fock/RPA theorjsi2i. We then 
use the full bosonized expression for the fermion propa- 
gator. Here we find large violations to Fermi liquid the- 
ory both at quantum criticality and in the nematic Fermi 
fluid phase. As an application we give a calculation of the 
fermion one-particle density of states. Finally, in section 
II VI we draw our conclusions. To help keep this paper self 



contained, in Appendix 1X1 we give a short review the ex- 
tension of bosonization to ZJ-dimensional Fermi systems. 
(For a more in depth review, see Ref.d). In AppendixlH 
we summarize details of the effective quadrupolar inter- 
actions, including fermion screening and Landau damp- 
ing effects. In Appendix O we discuss the effects of the 
(uncondensed) s-wave channel on the effective theory for 
the nematic. The details of the calculation of the boson 
propagators are given in Appendix IdI 



II. THE NEMATIC QUANTUM PHASE 
TRANSITION AND THE ORDER PARAMETER 



In this section, we consider the boson theory, obtained 
via an extension of bosonization to greater than one di- 
mension, near a nematic (Pomeranchuk) instability of a 
translationally invariant fermion system. In the nota- 
tion of Appendix^ we take the following action for the 
bosons 



^0 = 



N{0) 



/ (fxdt 



(2.1) 



and forward scattering interactions 




Here, S labels the patch defined by coarse graining the 
Fermi surface, and the density of quasiparticles in a patch 
may be obtained from the boson field (ps{'x.,t) via the 
relation 6nsix,t) = iV(0)vs • V.^5(x,i). Fs-t(x - x') 
is therefore the interaction between particle-hole pairs in 
patches S and T. 

We begin by analyzing our bosonized theory for a con- 
stant field configuration and reproduce Pomeranchuk's 
result. Consider configurations such that 6ng is constant 
in space and time over some particular range of time T. 
The resulting action is: 



VT 2 

5'shapo = "TTTTTY (1 + -^o) i^o) 



where we have expanded Sns as: 



N{0) 

VT 
2N{0) 



N/2 



5](l + F,)((m+)V(m7)') (2.3) 



'^"^ " V iv ^ [""^^ ^^^^ 

1=0 



9s) + mj sin{e9s)] (2.4) 



£>0 



and introduced the Fermi liquid parameters via 

Fs_t = ^Fo + ^^F,cos^(0s-^t) (2.5) 

e>o 

Hence, for arbitrary mf, we find that any Fi < —1 will 
destabilize the Fermi liquid. The point Fg = —1 we 
shall call the Pomeranchuk (nematic for £ — 2) quantum 
critical point (QCP). Though Fermi liquid theory breaks 
down at this point, Luttinger's theorem is still obeyed. 

It should also be noted that in the above analysis 
we could have included interactions involving large an- 
gle scattering (which may lead to charge density / spin 
density wave instabilities), corrections to the linearized 
dispersion, three- or four-body interactions and BCS pro- 
cesses in the bosonized theory. However, except for BCS 
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processes, these effects are irrelevant in the Fermi liq- 
uid phase though some become important near the the 
nematic critical point to be discussed below. 

For simplicity, in the rest of this section, we shall spe- 
cialize to the i = 2 instability in two-spatial dimensions, 
the 2D quantum nematic liquid crystal, though the re- 
sults generalize easily. 



the type 



^ +N/2 



ilBs 



(2.9) 



where ^ for ^ = {±2, ±6, ±10, . . .}. In terms of the 
mi's, the equation of motion becomes 



A. Saddle point expansion near the nematic QCP 

To study the nematic QCP, originally considered by 
Oganesyan, Kivelson and Fradkin^S, we set all to zero 
for ^ 2. This is reasonable since these other modes 
are not critical and their effect is only to introduce finite 
renormalizations of the parameters in the effective theory 
of the critical (quadrupolar) modes (see below). 

On the broken symmetry side (i. e. F2 < —1), the 
quadratic action is no longer stable. So, in order to make 
the theory consistent, we need (at least) a quartic term 
in the bosonized action. Here, as a specific example, we 
consider the quartic interaction that arises from correc- 
tions to the linearized dispersion in the bosonized form 
given in Ref.|36j, 



7iV(0) 



C 



(2.6) 



This term can be found by a direct extension of our 
e„ expansion of (|A12|I to third order, noticing that it 
is an exponential series in iSnsen/ {N {0)hvp) . An esti- 
mate of 7 may therefore naturally be obtained from the 
equations of motion of cither the boson or fcrmion pic- 
tures. However, this is not the only quartic contribution 
to the action since an 8-fermion interaction term (that 
is quartic in densities) also contributes at this level with 
a more general form similar to the Fermi liquid interac- 
tions. Nevertheless, this will give us a flavor of the broken 
symmetry phase. 

Since this effective theory is no longer quadratic in the 
bosonic fields, we will examine its behavior within a semi- 
classical approximation, which means that we will first 
find the extremal configuration and then expand about 
it. Thus, we write 



5ns = N{0)vs ■ 



(2.7) 



where Srig is the solution of the classical equations of 
motion at a uniform, mean field level: 



-F2(0)^COS 
T 



(,5ng)^=0 (2.. 



3!iV(0)^ 



where as stated earlier, we set F^(0) = for £ 7^ 2. Be- 
cause the non-linear term is cubic, we seek solutions of 



mi + F2(m2<5<?,2 + 'ti_2<5£,-2)+ 
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3!iV(0)2Ar 



mi^mi^me-e^-i^ = (2.10) 



£1/2 



We first observe that for F2 > —1, the only possi- 
ble solution is nii — for all £, as expected for the 
isotropic case. If F2 < —1, there exits a whole set of 
non trivial solutions involving, in general, all harmonics 
£ = {±2, ±6, ±10, . . .} obeying particle hole symmetry. 
Nevertheless, when we are near the phase transition, i. 
6- -^2^—1, one can find the set of solutions analytically 
with: 

m2m-2 = ^ [{m+y + (m^)') = -N{OfN \ l+F^iO) \ 

^ . (2.11) 
using the notation of (|2.4|l and for the higher harmonics 



|m4fe+2 |VA^(0)2iV -|1 + F2(0)| 
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2k+l 



(2.12) 



so that in the limit F2 — > — 1 we can neglect the higher 
harmonics. 

From this calculation, we conclude that near the 
F2 = —1 nematic QCP, the Fermi surface takes the sim- 
ple shape: 



7V(0)i 



^ 2|l + f2(0)| 

7 



cos(26'5 - a) (2.13) 



where a picks out the major axis of the ellipse that is 
spontaneously chosen. Without loss of generality, from 
now on we will set a = 0. We also conclude that fur- 
ther away from the critical point, the anharmonic quar- 
tic terms generically introduce higher harmonics to this 
shape. From this point of view for example, the Fermi 
surface in the nematic phase may become increasingly 
flatter away from the critical point leading to an addi- 
tional instability towards a smectic phase breaking trans- 
lational order in one direction. However, this phase tran- 
sition cannot be seen within our Forward scattering only 
model and all that happens deeper into the phase here is 
that the shape becomes more anharmonic. It should also 
be noted that the other harmonics, £ — {0,±4,±8, . . .}, 
appear when particle hole symmetry is broken, e. g. 
when a cubic term in the action is introduced (53) cor- 
responding to the addition of quadratic terms to the 
fermion dispersion. 
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B. Theory of the Quadrupole Moment Density 



and then switch over to the angular momentum basis: 



The next step beyond mean field theory is to formulate 
an order parameter theory. To that end, in this section 
we are interested in focusing on describing the behavior 
of the quadrupole moment density, mj(q, w). Again, for 
our present purposes we shall keep only ^2(9)- However, 
in Appendix ICl we show that by keeping also FQ{q), the 
additional effects of the non-critical modes do not change 
our results in any essential way. 

Now, our goal here is to obtain an action entirely in 
terms oimt ■ This can be easily accomplished by follow- 
ing Ref.jl6j and using a Hubbard- Stratonovich transfor- 
mation to aid the diagonalization. In the end, however, 
both the auxiliary fields and mf for £ 2 shall be inte- 
grated out. 

In momentum space, our free action can be written as 



5n 



(Pqduj 



-io] 



where 



X^(q,^) = ^(0) 



vs • q 



oj - 



vs ■ q 



= 7V(0)- 



cos 9 s 
— cos 6 s 



(2.14) 



(2.15) 



is the density-density response function in the "small q" 
limit of patch S with s ~ to/qvp. (Note: when needed, we 
regularize the denominator by letting s s + zesign (s) 
according to the usual time ordering prescription.) 

The interactions, described by the quadratic action 
•^int, c.f. Ea. (|2.2|l . become diagonal in the angular mo- 
mentum basis, i.e. in terms of the multipole densities 
mf{q,U!). In particular, the contribution to ^int from 
the £ = 2 (quadrupolar) densities, is 



'S'int — — 



qduj 



(27r)3 



/2(g) 



I riq{q,uj) I + I 7712 iq,Lj) \ 

(2.16) 

with f2{q) defined as usual in Fermi liquid theory through 
^2(9) = N{0)f2{q). The contributions from the other an- 
gular momentum channels have a similar form (in terms 
of the respective Landau parameters.) 

To aid the diagonalization, we split up the free part of 
the action, Sq, using the Hubbard-Stratonovich transfor- 
mation: 



1 f cPgdio 



Xsil> w)cts((?, uj)as{-q, ~uj)~ 



(as{q,uj)Sns{-q,-oj) + (Ts{-q, -u})Sns{q,uj)) 



(2.17) 



Sn = - 



d^ qdbj 
~{2^ 



In the large- limit, we have 



N{0) 



2tt 



dO cos 6 
2ns — cos 



(2.18) 



(2.19) 



Here we have Fourier transformed with respect to q = 
(cos 0, sine/)), that is, 9 = 9s — (f). 

Now, we will integrate out all the densities, ex- 
cept for £ = 2. This can be easily done since 5*0 is a 
linear function of these fields, while they are absent in 
S'int. This is so for this model with only a quadrupolar 
interaction , i.e. we have set their corresponding Fermi 
liquid parameters to zero. (In the vicinity of the nematic 
transition it is straightforward to include the effects of 
the I ^ {) channels. Their net effect is to give rise to sim- 
ple renormalizations of the effective theory we are about 
to derive. A detailed analysis is given AppendixEJ. The 
result is a delta function for the ag fields, allowing us to 
also integrate them out with the net result to simply set 
cr^ = for all I except t — 2. This gives us the following 
simple expression for the effective free action: 



1 /" d'^qduj 

^o=-2^y 72^ 



ix'o + Vxl) 



2 (a. 



(2.20) 

and we now have an action entirely in terms of the 
quadrupole moment density. 

The final step is to integrate out the af fields. This is 
easily accomplished and we obtain the Gaussian level of 
the order parameter theory, including the effects of the 
interactions S'int. The action of the effective theory is: 



2 J [^T^r 1X2 X2 

(2.21) 

where X2 (*i 1) dynamical correlation function (sus- 

ceptibility) of the quadrupolar densities: 



x^(s,g) = 



x'ois)±xlis) 



l-f2iq){x"oi^)±x'M 



and 



X2i = N{0) 



- Se,o + Ko{s) 



1 - Ko{.s) 
l + Ko{s) 



(2.22) 



(2.23) 
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with 



Ko{s) 



(2.24) 



Naturally, this is just RPA quadrupolar susceptibility of 
Ref. 

We should stress that the effective action S2 of 
Ea. (|2.21|l does not include the effects of the non-linear 
interactions represented in 6*4, c.f. Ea. H2.6() . which are 
crucial to stabilize the nematic state past the nematic 
QCP (for ^2(0) < —1). As we discussed above, these 
non-linear terms do mix the different angular momentum 
channels. However, provided there is no condensation for 
£ 7^ 2 it is still possible to integrate out these degrees of 
freedom, at least perturbatively. Thus, sufficiently close 
to the nematic QCP, the effects of the higher angular 
momentum channels will remain perturbatively small. 

Furthermore, should we have decided to integrate out 
the density fields, m J in favor of the auxiliary fields ti^ , 
we would have found the propagators of the af fields to 
be the RPA effective interaction: 



hiq) 



l-/2(9)(xS](s)±x5(s)) 



(2.25) 



The action of the fields is precisely that found by 
Oganesyan and coworkersi^ also obtained through a 
Hubbard- Stratonovich transformation but performed di- 
rectly in the fermion theory. We shall find that under- 
standing this interaction is the key to understanding the 
physics of the nematic QCP. 

C. Order parameter theory of the nematic QCP 

The theory with the action given by Eas. H2.f| ) and (|2.2|) 
seemingly describes the quantum mechanics of a fluctu- 
ating surface. This suggests that the effective degrees of 
freedom ought to be long-lived bosonic modes of the fluc- 
tuations of the shape of the Fermi surface. In other terms, 
this bosonizcd theory would seem to be entirely described 
by stable collective modes. However, these bosonic ex- 
citations are not generally stable due to Landau damp- 
ing effects, represented by the branch cut singularities in 

19.41 



We will examine this problem more closely. It is useful 
to introduce the quadrupole density spectral functions: 



S2{q,s) 



where B2 {q, s) denotes the spectral function of V.^ {q, s), 
the correlation function of the fields. 

The analysis is greatly simplified upon recognizing that 
V^{q,s) is a polynomial function of K^){s). For £ = 2 



-2 sign {s)Imx2 {q, s) 
2 sign {s)JmV^ {q, s) 



B}{q,s) (2.26) 
f2{qy 



S+iq,s) 



FIG. 1: The (normalized) quadrupole density spectral func- 
tion S^{q,s) = S:^{q,s)/N{0), Eq. ijOfit and Ea. l^^ for 
g 7^ 0, very close to the nematic quantum phase transition 
from the Fermi liquid phase, -^2(0) — > — 1+. As the QCP is 
approached, there is a large increase of spectral weight in the 
the 2 = 3 overdamped quadrupolar mode at very low fre- 
quency. Notice the delta function contribution of the propa- 
gating mode with z = 1 discussed in the text. The propagat- 
ing mode with z = 2 gives a similar delta function contribu- 
tion (with smaller spectral weight) to ^^(g,^). 



Elwe show that V^{q, s) have the partial fraction expan- 



V2^iq,s) 



(q) 



(2.27) 



We 



where (3 — a,b,c for ¥2^, and /3 — a,b for V2 
can therefore view V^{q,s) as a sum of terms each of 
the form of the RPA s-wave channel effective interaction, 
renormalized by a residue and with an effective inter- 



action {q) 



(1-^^(1))- 



Details of these expansions 

are given in Appendix iBl 

Near the nematic QCP, where 1^2(0) 
defining the quantity 52{q) 



I| < f, upon 



S2{q) 



1 + j^2(g) 

F2{q) 



= 1 



1 



i^2(0) 



Kq 



(2.28) 



the effective interactions V^^ (q, s) become simple and we 
obtain: 



V2+{q,s) 
V2-{q,s) 



1 



2N{0) 
1 



1 



1/4 



Mi) 
2 

f 



Kois) 



1/2 



47V(0) , 



Mi) 
4 



(2.29) 
(2.30) 



Thus, near the critical point we find a propagating mode 
with dispersion uiq — (l/^/2)qvF {b, z ~ 1 mode), 
represented by the pole in the second term for 
(see Ea. (|2.29ll l. and another one with dispersion ojq — 

(^V^^iq)/^^ qvF (a 2; = 2 mode), given by the pole in 



the denominator of 1/^2 (g, s) is a cubic function of i^o, ^2 (^^^ ^"-l- 



Furthermore, we also find an over- 



whereas for V2 {q, s) it is quadratic in Kq. In Appendix damped mode, given by the pole in first term of (see 
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Ea. (|2.29|l '). with a dispersion relation of s = —i52{q)/2, 
and dynamic critical exponent z = 3. For q sufficiently 
small, and asymptotically close to the nematic quantum 
critical point, we find that this overdamped mode domi- 
nates the spectral function over the other two (propagat- 
ing) modes. 

In Fig. n we present a plot the spectral function 
i3^(g, s) = — 2 sign (s)XmV2^(g, s). It shows that, as the 
QCP is approached from the Fermi liquid side, there is a 
large transfer of spectral weight in the quadrupolar spec- 
tral function to the low frequency end of the spectrum, 
associated with the emergence of the overdamped 2 = 3 
mode. This mode thus controls the quantum critical be- 
havior. 

Conversely, normal Fermi-liquid behavior is obtained 
if i52(0) is finite so that z = 1. On the other hand, if the 
overdamped mode were to be absent, the dynamic crit- 
ical behavior would be controlled by propagating mode 
with z — 2 discussed above. The trend is thus opposite 
to what one might naively expect: the higher the z, the 
higher the effective dimension D + z, the stronger the di- 
vergence of the overdamped mode. Recently, Yan^^ pro- 
posed that the transition to the quantum nematic state 
should have dynamic critical exponent z = 2. The anal- 
ysis we just presented shows that this is not the case. 

As a result of the above analysis, we obtain the follow- 
ing effective action for the quadrupole density near the 
nematic QCP: 



'S'qcf 



1 



(Pqdijj 



{2i\s\-l~F2{q))\m+\ 



2iV(0) J (27r)3 

(4.2_i_F2(<z))|m-p] 
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d'^xdt(m^'' + m^'^y (2-31) 



8iV(0)3A^ 
The order parameter field is 



s 

'•"^(q^^) = \I^^Snsi(l,uj)sm2{9s - 



(2.32) 



where, again, cj) is the direction of q = q{cos(j),sm4>). 
This action is therefore very similar to Hertz's action for 
the ferromagnetic quantum phase transition in itinerant 
fermionic systems^-, but here within the context of the 
formation of nematic order (and similar actions may also 
be obtained for higher £). Recently, Nilsson and Cas- 
tro NetoiJ derived this action using Fermi liquid theory 
methods. 

Now, let us look on the broken symmetry side, in the 
nematic phase. As in an in any theory with an 0(2) 
symmetry, here we will find that the order parameter 
will spontaneously pick a direction. As a result, it is no 
longer useful to Fourier transform with respect to 4> in 
Eq. (|2.32|l . Rotating back and after the saddle point ex- 
pansion about the classical configuration H2.13|l . we find 
the quadratic action on the nematic side: 



^nematic 
'QCP 



d'^qdbJ 
2N{0) J (27r)3 



m2(q,w)-X2 ^■m2(-q, -w) 



(2.33) 



where m2(q, cj) = (toj (qi "^2 (qj^)) ^'^'^ 



x2' = m-^<i^^\^+F2m){^Q J) 

/cos4^ sin4j\ A 

' ' \ sm40 —cos 40/ ' ^ \ — 1 



(2.34) 



in the reference frame in which the nematic order param- 
eter is diagonal, i.e. its principal axes, whose orientation 
is determined spontaneously. The difference between this 
action and the previous one discussed above for the sym- 
metric phase, is the emergence of the last term which 
originates from the non-linear (quartic) term in the ef- 
fective action. This term ruins our ability to rotate (j) out 
of the action. 

Noting that the off-diagonal terms are higher order in 
s, we may write this in the simplified form: 



nematic 
QCP 



d'^q 

2N{Q) J (27r)2 



duj 
2^ 



(2i|s|cos2(2(^) - 2|l-t-F2(0)| - Kq^)\m+\^ + (2i|s| sin2(2(^) - Kq^)\m-^ 



(2.35) 



We see that now is the amplitude mode, which has on the broken symmetry side has: 



2 = 1, while TO2 is the nematic Goldstone mode which 



continues to have dynamic critical exponent z 
in the nematic phasei^ 



3 even 



The last point is to discuss how this affects the original 
boson theory. If we bring back all the integrated out 
angular momentum channels, we find that the free action 



X^(q,L.)=A^(0)- 



vs ■ q 

- vs • q 



with a weakly renormalized Fermi velocity 

(1 + 4| 1 + ^^2(0) I COS^ (205)) V5 



(2.36) 



(2.37) 



8 



Hence, spontaneous symmetry breaking essentially pro- 
duces a Goldstone mode that continues the critical, 2 = 3 
behavior into the nematic phase while leaving the rest of 
the theory virtually untouched until deep into the broken 
symmetry phase. 



III. FERMIONS IN THE CRITICAL REGIME 

In the past sections we discussed the behavior of the 
collective modes near the nematic quantum phase transi- 
tion. Much of what we discussed in the previous section 
on the behavior of the collective modes is indeed in com- 
plete agreement with the RPA treatment of this theoryi^. 
This should not be a surprise since RPA is asymptotically 
exact at low energies and at low frequencies. This is also 
the reason while bosonization works in the same regime. 

We will now turn our attention to the behavior of the 
fermionic degrees of freedom near the nematic QCP and 
and in the nematic phase. This is very different problem. 
In Ref.To'l the behavior of the fermion Green function 
was studied pcrturbatively and a startling non-Fermi liq- 
uid behavior was found already at the lowest ("Fock") 
order. However, this very finding raises questions on the 
applicability of perturbation theory for the fermion prop- 
agator. In this section we will use bosonization methods 
to address this problem. 

Here we will use bosonization to compute the fermion 



propagator. Within this approach one has a theory for 
the bosonized degrees of freedom and a set of operator 
identities relating observables of the fermionic theory to 
those of the bosonic theory. For a summary see Appendix 
IXI There are two important issues to keep in mind. One 
is that the bosonized theory is exact for a fermionic the- 
ory with a linearized dispersion and forward scattering in- 
teractions {i.e. those described by Landau parameters). 
The other is that one has, within this theory, an oper- 
ator to represent the fermion. Corrections to the linear 
dispersion as well as other (non-forward scattering) inter- 
actions are represented by non-linear terms in the action 
of the bosonized theory. The expressions that we will 
derive below apply strictly speaking to the fixed point 
theory, in which these perturbations are not included. 
It turns out that corrections due to the non-linearities 
of the fermion dispersion and other such terms are ir- 
relevant (both in the Landau phase and at the quantum 
critical point). As such they will affect the results at high 
energies and at momenta but their effects become negli- 
gible in the low energy limit. Please note that one such 
operator, discussed in the previous section, stabilizes the 
nematic phase, i.e. it is a prototypical dangerous irrele- 
vant operator. 

The boson Green function may be found from our ac- 
tion in a similar way as we found the above density- 
density correlation functions, that is, using a Hubbard- 
Stratonovich approach. The result is: 



J 



(5,T)(x,t)=GO,(5,T)(x,0+* j |^G°,(s)(fc,e)V^s,T(fc,£)G^(^)(fc,£)(e^C'--*)-l) (3.1) 



function: 



where G^^^^^ is infrared divergent unless x || and where the relative angle, 9s — 4> appears here (see Eq. 
S^T. On the same patch, G%^^s s)(^^*) given by the EEl)- We also have used in ^ the free fermion Green 
standard expression 

G?5(5,S)(X'0 = (¥'5(x,t)^5(O,0)) - ((^5(0, 0)2) 
' fis ■ X + ivpt + tasignt 



In 



la 



Gf(s)(9'^) 



LU — -vs ■ ci + ie sign (w) 



(3.5) 



(3.2) 



tor (see Eq. I)A20(I of Appendix to give us the fermion 
Green function of the form: 



where a is a short-distance cutoff. In Eq.JHJ we have We will now useth e bosonized of the fermion opera 
denoted by Vs.xil,^) the effective interaction 

^E^^^'''^"'^''''^'^"'"(-^(9'-)-H-9,--)) 



Gi.(x,<)-E (^) {T yjsi^,t)i;!,iO,0))e^--- (3.6) 



(3.3) 



^EGF(s)(x,t)e»''-'' 



which for the quadrupolar case, becomes simply 

VsMq^^) = V+{q,io) cos2(0s - 0) cos2(0T - </>) 
-I- V~ (q, Lu) sin 2(03 - (f>) sin 2{9t - (f) 



N 



for which we find the explicit expression 



(3.7) 



(3.4) 
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G'F(s)(x,i) 



cPkde 



,i(k-x-et) 



(3. 



r 



This expression has many similarities with the bosoniza- 
tion formulas usually obtained in one dimension. In 
particular, the free-fermion pre-factor also arises there. 
However, the exponential factor, which in one dimension 
yields an anomalous dimension for the fermion operator, 
plays a very different role in dimensions higher than one. 



Ref.|3- To do that we will expand the exponential and 
Fourier transform to momentum space to find its dia- 
grammatic expansion. The result will be a series of con- 
volutions since in real space they are products. The first 
order term is: 



A. Diagrammatic expansion for the bosonized 
theory 

We will first show that the bosonized formula of 
Ea. (|3.8|) is consistent with the perturbative results of 



J 



(27r)3 



F(S) 



(q-k. 



VsMk,e) 



(3.9) 



which does not look like it obeys the Feynman rules 
for the perturbation theory of non-relativistic fermions. 
However, we may utilize the following identity 

(G°,(5) (q - k, - e) - G°,(5) (q, CO)) G%^s) (k, s) 

= G°,(5)(q-k,w-e)G^(5)(q,c.) (3.10) 

and use it again through 
G'F(s)(q-k,t^-e) = 



G%is) (q - k, w - e) - G^(s) (q, uj) + G"p^s) (q. ^) 



(3.11) 



To obtain 



SG^P{q,Lo) 



G%s)i'l,i^) 



(Pkde 

^ ( G^s) (q - k, - e) + G^p^s) (k, e) ) ^5,5 (fc, e) 



(27r)3 



(3.12) 



The second term is actually the shift in the chemical 
potential, S(A:^,0), and is zero by the effective particle- 
hole symmetry of this theory (with a linearized fermion 



dispersion). We therefore obtain, 



SG^P{q,io) 



GF{S)il^^) 



(fkde 



G°,(s)(q-k,L^-£)l^s,s(fc,e) (3.13) 



which is the correct result since, as usual, the Hartree 
term vanishes. 

Using the same tricks and assumptions, the second or- 
der contribution can also be worked out (though it is 
much more work). The series to second order is shown 
in Fig. 121 please keep in mind that the interaction, Vs.s 
is the full bubble summed RPA interaction. Higher di- 
mensional bosonization therefore keeps all diagrams in 
perturbation theory that contain up to simple bubbles 
while neglecting more complicated bubbles as is usual 
in the RPA (in the Landau theory of the Fermi liquid, 
these irrelevant operators contribute subdominant poten- 
tially non-analytic temperature and frequency dependent 
terms to physical quantities'^). 

Some time ago, Kopietz and Castillck^SiS used a some- 
what different (and in principle equivalent) form of 
bosonization and discussed the effects of a quadratic term 
in the fermion energy dispersion. However, instead of 
using an operator identity (and thus not exploiting the 
non-perturbative character of bosonization) they chose to 
make contact with perturbation theory, and proceeded to 
propose a modified form for the fermion propagator di- 
rectly. It gave an exponential factor similar to the one 
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theory 




FIG. 2: Bosonization's Feynman diagram series. 



appearing here except the fcrmion Green functions that 
appear in the exponential include the quadratic terms in 
ther energy dispersion. However, these Green functions 
don't satisfy the criterion E(A:f,0) = so that in ad- 
dition to the exponential, they include a pre-exponential 
factor that is necessary, for example, to cancel the second 
term of eq. (|3.12() . As a result, order- by-order in Vs^s, 
one needs to keep precisely the right pre-exponential fac- 
tor to cancel the additional terms. A calculation with 
their method can therefore only be carried out to a small 
finite order, and one can no longer think of the exponen- 
tial factor as separate from the pre-exponential factor. 
In contrast, we have seen here that the bosonized expres- 
sions, when treated consistently, yield exact results which 
agree with those of perturbation theory order by order, 
albeit only in the low energy limit, including the singular 
behavior. 



B. Perturbative results 

Before computing the full non-perturbative form of the 
fermion Green function, let us verify that our bosonized 
theory, reproduces the perturbative results of Oganesyan, 
Kivelson and Fradkin^^ near the nematic QCP. We shall 
be interested, therefore, in the integral: 



(3.14) 

Here Vss{k, e) is the effective interaction mediated by the 
collective modes, c.f. Ea. lfT^ . Ea. lfOHl and Ea. l|njl . 
The main contribution to the self-energy is due to the 
overdamped Goldstone modes, as noted in Ref. 19]. Thus 
we take a generic interaction of the form 



Vs,sik,e) = 



1/N{0) 



6ik) - Kois) 



(3.15) 



noting that cos^{2{9s — 0)), which should appear as a 
coefficient in the nematic case, only introduces irrelevant 
contributions to the integral and is therefore left out of 
this analysis. Here, at the QCP, d{k) = —Kk'^/2 but as 
a check we may characterize the Fermi liquid phase by 
letting 6{k) = S, a constant. 

Before computing S]g''(q, w), we should note that this 
expression gives the clearest definition of the patch. 
Given the ultraviolet cutoff kp — A/2 < |k| < kp + A/2, 
the patch is defined so that in comparison to the original 



I](i)(k = ks + q,c^)«4'^(q,c^), qeVs (3.16) 



and therefore the patch width is of order 



A V^!fA 



(3.17) 



dictated by the curvature of the Fermi surface alluded to 
earlier. 

Now, to simplify our calculation, let us focus on the 
fermion lifetime near the Fermi surface: 



2Tsiq,iu) = -2sign(tj)XmEs(q,t^) 



(3.18) 



This can most easily be expressed in terms of the spectral 
function, _B2(q, s), derived from the imaginary part of 
Vs.siq, s) as in Ea. lfT^ and Ea. (|0) (after setting the 
angular factors to 1) 



2rs(q,c^) 



-A/2 27r Jq,^-\^\/vp 



dkn ^ A UJ - Vpign - kn) 

2tt \ ' kvp 



(3.19) 

where qn='^s ■ i/^^f- Here we note that this integral is 
dominated by the contribution of the overdamped mode, 
which enters in -B^(q, s). (The other contributions, asso- 
ciated with the propagating collective modes, only yield 
regular dependences in the fermion frequency.) Setting 
Qn — > and looking at the limit a; — > we obtain 

Tsico)^-^ dk^ dhB+{\lkl + klf) 

(3.20) 

using normal and tangential coordinates. From here, we 
first integrate over the tangential momenta, look at small 
q„ and obtain the long wavelength behavior. We find that 
in the Fermi Liquid phase {5{k) — const) 



rs(c^) 



1 f UJ 



— I \oguj/KvF\ 



N{Q) \VF, 

while at the nematic QCP, using 5 = — Kfc^/2: 



.(- 



N{Q) \kvf 



2/3 



(3.21) 



(3.22) 



In the Fermi liquid, the energy of the quasiparticle is pro- 
portional to Lu and we classify them as long lived. How- 
ever, at the nematic QCP, one can show via Kramers- 
Kronig that the real part of the self energy also goes like 
jcjp/'^. Hence at the nematic QCP, the life time of a quasi 
particle is not well defined and to try to understand it 
as a perturbation about a free theory of long lived quasi 
particles is meaningless. 
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Thus, the bosonized theory reproduces the results 
found earher on perturbatively by Oganesyan and 
coworker^iS (see also R,ef . poEoj . It should be noted 
that the uP'l'^ law was also found to appear in the pertur- 
bative calculation of the fermion self-energy in a model 
of holes interacting via a forward scattering C/(l) gauge 
interaction with a similar form to l|3.15|l in the context of 
high temperature superconductors^^-^'^, and in the per- 
turbative treatment of the quantum critical point in a 
ferromagnetic metali^ 



C. Non-perturbative results 

We now return to the non-perturbative bosonized ex- 
pression for the fermion propagator of Eq. (|3.8() , which we 
will write as 



G 



.^^s)M ^ Zs{^.t)G\s)M (3.23) 

In the Fermi liquid phase and at long distances and low 
frequencies the factor Zs approaches a constant value, 
Zs — Zp < 1, i.e. the quasiparticle residue of the Fermi 
liquid state. Our goal here is to investigate the behav- 
ior of Zg{'K, t) near the nematic QCP and in the nematic 
phase. However, given the complexity of the full analytic 
expression, in this paper we will consider only on the 
equal-time, t = 0, behavior (sometimes called the "one- 
particle density matrix") and the equal-position, x = 0, 
dynamical correlation function, and only at zero temper- 
ature. We will discuss its full behavior elsewhere. 



The 



time fermion propagator 



As in the perturbative calculation, it is convenient to 
express In Z as an integral over the spectral functions 
Bf(q,s). Once again, at the nematic QCP, the im- 
portant contribution is due to the overdamped collec- 
tive mode in -B^(q, s), and we will neglect all other con- 
tributions. (This is an approximation which gives the 
long time behavior accurately. An expression valid for 
all times also includes the contribution of B2 {q, s).) 

Since we have set i = here, the result is quite simple 



lnZ(x„,0) = / -^/2(fc„)(cos(fc„a;„) - l) (3.24) 



where 



A/2 27r Jo 27r (cj + ufA:„)' 



We find, both from a numerical computation and from 
an analytic estimate, that 

1 (I. \ 1 h {knVK) „„^ 



where f2{kn\fi^) is a regular function of its argument. 
From this analysis, we conclude after performing the final 
Fourier transformation: 



Zs{xn,'^) =Cpexp. 



1/3- 



(3.27) 



vahd for |a;„| » ^/n. Here b — 0.0658 and Cp is 
a constant factor resulting from subdominant terms in 
hi^n)- This sharp decay of Zs{Xn,0), faster than any 
power law, introduces a scale (similar to a correlation 
length arising from a gap in the spectrum) and dom- 
inates over the Fermi liquid behavior at low energies. 
From the above expression, the correlation length is of 
order ^ ^ ^/K,(kp^/K)'^ which is much longer than length 
of the interactions, ^/k. 

Let us also compute the behavior of Zs on the nematic 
phase by focusing on the effect of the Goldstonc modes. 
Recalling our discussion of the order parameter theory, 
we replace ()3.15|l with 



VsAk,e) 



sin2(20s) 



1 



A^(0) -Kk^/2 + i\s\sin^{2<j)) 



(3.28) 



which is the contribution from the Goldstone mode 
{Kq{s) k, —i\s\). In the limit, fc„ 0, that is on the 
Fermi surface, ^5 4- 7r/2 and therefore, we simply 

have 



Vs,s{k,e) 



1 



N{Q) V sin2(20s; 



(3.29) 



so that the difference between this and the symmetric 
side is simply that n k/ siv?{29s). Hence, we may 
directly write down Zs'- 



Zs{xn,^) = Cp{^s) exp' 



6|sin(20s)| 



4/3 



1/3- 



N{Q)vF^ 

(3£0) 

where h is the same constant of Eq. (|3.27|) . In Eq. I|3.30(l 
we have not included the subdominant contributions 
which become the leading terms along the symmetry- 
dictated directions, the "nematic axes" , along which the 
angular factor vanishes. Along the nematic axes the be- 
havior of the equal-time correlation function has a more 
Fermi liquid like long distance behavior as shown by 
Oganesyan et. al., but due to the introduction of patchs, 
we cannot accurately capture this behavior here. 

Thus, we see that on the broken symmetry side, we 
have special points at ^5 — n7r/2 where the Goldstone 
mode weakens and subdominant behavior takes over. At 
these points, Zs — Cp(^s) — < 1 and the quasiparti- 
cles become long lived. It is interesting to note that the 
angular dependence shown here with a power of 4/3 is 
similar to that of the perturbative calculation using the 
same transformation of k ^ k/ sm^jWp) on Eq. 
and it agrees with the results of Ref . |l3 • 
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2. The fermion residue 

One simple calculation we can do with the above result 
for Zs{x„,0) is the fermion residue following Migdal^. 



Zg = n(ki? - q) - n(kF + q) 



(3.31) 



(3.32) 

where the exponential factor tells us to close the contour 
in the upper half plane. This expression may be written 
in terms of the real space, real time fermion Green func- 
tion, at time t = — e. Inserting our expression for Gp 
within Bosonization and neglecting interpatch scatter- 
ing (which should produce analytic in q/A contributions 
here), we obtain: 



Z„ 



smqnXn 

dXn Zs (X„ , 0) 



(3.33) 



In the Fermi liquid phase, we find Zs{xn,0) = Z < 1 in 
the long distance limit and this leads directly to Zg^ — Z 
when g„ — > as expected. 

At the nematic QCP and into the nematic phase away 
from the nodal points discussed earlier, Zs{xn, 0) is short 
ranged and we may expand to leading order in g„ and 
obtain 



12 . 

— 9n4 
TT 



(3.34) 



with the correlation length ^ — ^/^{^/KkF /'iirhY . Thus 
the fermion residue vanishes linearly similar to how it 
would in a Fermi liquid at finite temperature where tem- 
perature makes the correlations short ranged. It should 
be noted, however, that because Zs(a;„,0) decays slower 
then e^"^" in the long distance limit, this series expan- 
sion that we have used is poorly defined at higher order 
with the coefficient of q^'^'^ growing so rapidly that the 
Taylor expansion has zero radius of convergence in the 
complex-(7„ plane. Thus Zg^ is not analytic in q„ and 
the Fermi surface may still be defined as a singular point 
in n(k). The same behavior occurs in the nematic phase 
for generic momenta, except along the directions of the 
nematic principal axes where a finite residue is obtained. 



In contrast with our result for the equal-time correlation 
function, Ea. (|3.35|l approaches a constant at long times. 
However, it decays to that constant much more slowly 
than in a Fermi liquid where we would expect the expo- 
nent 2/3 to become 2. As a result, the non-perturbative 
effects are here less important and, consequently, this 
time-dependence appears to exhibit the same power-law 
behavior as the life time calculated perturbatively in Sec. 
ElBl 

The coefficient of the exponential was found to be 
A oc l/(fcp'AK) and deserves some attention. In a Fermi 
liquid, we would find A oc X/kp <C 1 and the time 
dependence would present a small correction from the 
Fermi-liquid behavior. However, A need not be small if 
1/k ^ kpX ~ A^, where A is the patch width (c. /. Eq. 
(|3.17|) '). Due to the emergence of A here, this limit oc- 
curs precisely where the Fermi-surface curvature begins 
to matter, and where the interaction length scale is still 
quite smah {kp'^ < < A"^). 

In the nematic phase, again letting k k/ sir? 29s, 
we find 



Cp(5) exp \ -A{\-iVi) |sin26ls 



4/3 



In I Sin 20s 



[vpt/^f'' 



(3.36) 



so that the angular dependence is similar to that of the 
equal-time behavior. 



4- The one-particle density of states 

We close this section with an application of these re- 
sults to the calculation of the one-particle density of 
states (DOS) in both the Fermi liquid and the nematic 
phases, and at the nematic (Pomeranchuk) quantum crit- 
ical point. 

The one-particle DOS is defined by the standard ex- 
pression 

7V*(w) = -sign(cj) - Im Gf{x,x;oj) (3.37) 



3. The fermion auto- correlation function 

This case turns out to be more complicated than the 
equal-time expression and we present a full analysis in 
Appendix^ By a simple integration-by-parts, we found 
that the double pole of Eq. (|3.25(l may be reduced to 
a single pole and that the integrals involved were less 

1/3 

singular by a full power (diverging like 1 / fc„ unlike Eq. 
(|3.26l) ') but with logarithmic corrections. From that anal- 
ysis, we found the general form of Zs to be 



Cpexp|-A(l-^V3) ^"^""^/^U (3.35) 



In Fermi liquid phase, using Zs(0, i) = Cp, a constant 
dependent upon the Fermi liquid parameters that goes 
to 1 for the non-interacting case, as expected we find 

Gp(x,x,L.) ^^Icp / dtG%s){0,t)e'^' 



s 

^Cf I dt 



N{0)vF 



vpt -\- zasign (t) 
= -in sign {uj)CfN{0) (3.38) 

so that here, N*{uj) — CfN{0) as expected. Note: Cf 
is not the fermion residue and may be greater than one. 
At the nematic QCP, we have instead 
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F(S) 



(0,i) 



(3.39) 



r 



so that N*{lu) = CpN{0)I{u) with 



-Tie 

./O 



■ sm (wu j exp 



-A(l-i\/3) 



(3.40) 

with a) = ^/kuj jvp^ vahd for uj <^ vp min{K~^/^, A}. No- 
tice that we have used throughout these expressions only 
the long time limit of the exponential factor. At shorter 
times the behavior of the exponential should be dom- 
inated by high energy effects which are insensitive to 
whether the system is in a Fermi liquid, a quantum criti- 
cal point or in a nematic phase. Thus, the time integrals 
have an implicit short distance cutoff (which we have de- 
noted by "0" ). In any event, we are only interested in the 
low frequency behavior which is dominated by the long 
time part of the integration range. 

By inspection we see that the exponential factor ap- 
proaches unity (quite rapidly) for large u 3> 1. Thus, 
the main effect of this factor is a correction to away 
from its value at zero frequency, i.e., /(O) = 1. The lead- 
ing finite frequency behavior, as w 0, is obtained by 
expanding the exponential factor in Eq. 



1 



2A 



du sin uju 

2/3 



Inu 



u u 



'KU! 



Vp 



2/3 



In 



Vp 



(3.41) 



where the last line is accurate for y/KUj/vp ^0.1. The el- 
lipsis in Ea. H3.41|l represents subdominant contributions 
at low frequencies, which vanish faster than uj'^/^lnuj as 
Lu ^ 0. As a result, the w^/'^ behavior of the inverse life 
time T^uj) calculated perturbatively, appears here as a 
cusp in the DOS (with a logarithmic correction). Unlike 
the lifetime, though, here A depends on the product of 
the patch width cutoff y/kpX ~ A and ^/k while r(w) 
depends only on the product of kp and ^/k. 

Thus we find that the low frequency one-particle den- 
sity of states has the form 



N*{uj)^ N* (0) + B In + 



(3.42) 



where 7V*(0) = CpiV(O) > iV( 0), si nce we found 
the constant Cp > 1 (see Ea. (1017^ 1. and B = 

ACpA^(0)^r(l/3) (see Ea. dDTSt 'l. Hence, the zero fre- 
quency value of the density of states is larger than the 



Fermi liquid value. As the frequency increases, the one- 
particle density of states decreases from its zero frequency 
value according to the uP'l'^ Inw correction term. This is 
a cusp singularity at a; = 0. It is important to stress 
that we obtained these expressions upon expanding the 
exponential factor in the auto-correlation function. This 
is consistent since this factor asymptotically (and rather 
rapidly) approaches 1 at long times. This also implies 
that, in this regime, our results should be consistent with 
the behavior of the fermion Green function found in per- 
turbation theory around IIartree-Fock/RPAi2i22i2Ma. 



In the. nematic phase, an angular average of Eq. 
enters our expression for N*(uj). In the long time limit, 
when the argument of the exponential is much less than 
one, we expect that for low frequencies 



A 



sin 263 



4/3 



A w 0.58A 



(3.43) 



Hence, in the nematic phase, the Goldstone modes con- 
tinue the critical behavior of this function in only a mildly 
weaker form. 

In summary, in this section we used bosonization to 
compute the non-perturbative behavior of the fermion 
propagator using the bosonized form of the fermion op- 
erator. We first checked that the non-Fermi liquid behav- 
ior of the nematic QCP and in the nematic phase, which 
were obtained earlier using conventional diagrammatic 
(perturbative) methods^Si^ is recovered here upon ex- 
panding the bosonized expression to leading order in Vss, 
i.e. a single boson exchange. However, upon a closer 
examination of the full bosonized result we found that 
the equal-time fermion correlation function has a much 
more singular behavior that could have been predicted 
in perturbation theory. In contrast, the fermion auto- 
correlation function (and hence the one -particle density 
of states) is seemingly consistent with the perturbative 
analysis of the quantum critical behavior. 

We note here that Chubuko\iS, has analyzed the quan- 
tum critical behavior of ferromagnetic Fermi liquids and 
claims that the uj'^^^ behavior, found at lowest order, per- 
sists to all orders in perturbation theory. The results of 
this section for the fermion auto-correlation function ap- 
pear to agree with those of Chubukov and coworkers. 
However our results for the equal-time fermion corre- 
lator apparently disagree with these results. Clearly, a 
more detailed analysis of the bosonized expression for 
this propagator is warranted. We will discuss this prob- 
lem in a separate publication. 
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IV. CONCLUSION 

In this paper, we have utilized the method of high di- 
mensional bosonization to study non-perturbatively the 
quantum phase transition from a Landau Fermi liquid 
state to a nematic phase, a nematic (Pomeranchuk) in- 
stability. For this purpose, we have constructed an order 
parameter theory from the boson theory by integrating 
out non-critical modes and verified that this boson the- 
ory is equivalent to RPA. We then turned to studying the 
bosonization form of the fermion propagator and found 
its diagrammatic expansion proving the correctness and 
clarifying the arguments leading up to that expression. 
This diagrammatic expansion keeps all diagrams up to 
the simple bubble in the spirit of RPA as applied to 
the density-density propagator and shows, in specific, 
that bosonization goes beyond the self-consistent Born 
approximation to include vertex corrections. We then 
found explicitly that bosonization reproduces the results 
of Hartree-Fock with an RPA interaction by showing that 
the lifetime computed in this limit also has an jwp/'^ de- 
pendence as originally found for the case of the nematic 
QCP by Oganesyan and coworkersiS. 

Lastly, we calculated the fermion propagator non- 
perturbatively and found the dramatic effect of the 
overdamped critical mode that induces short ranged 
spatial correlations that decay nearly exponentially 
^g-const. \x\ while the auto-correlation function ex- 
hibits a milder non-Fermi liquid behavior of the form 
jij- exp(— const. From this short ranged be- 

havior of the equal time Green function, we verify that 
the fermion residue vanishes at the critical point and 
into the nematic phase except at four special points. 
We also calculated the one-particle (fermion) density of 
states N*{uj). We found the low frequency behavior 
N*{uj) = N*{0) + B uj^/^lnuj (with N*{0) > N{0)) both 
at the quantum critical point and into the nematic phase. 
Thus, the fermion propagator exhibits unexpected be- 
haviors which could not have been anticipated by the 
existing perturbative resultsi^iiSiSfliMiSi. In a separate 
publication we will present a more detailed analysis of 
the fermion spectral function in both phases and at fi- 
nite temperature. 

Two recent papers, one by Yang'^* and another by Nils- 
son and Castro Neto^, also study Pomeranchuk nematic 
instabilities in two-dimensional Fermi systems. Yang also 
derives an order parameter theory within high dimen- 
sional bosonization. However, contrary to our results, 
concludes that the critical mode has z = 2 and is un- 
damped. While we agree that a z = 2 propagating mode 
does exist at the critical point, we find that the spec- 
tral function is completely dominated by the overdamped 
z — 3 mode. This effect is due to Landau damping, a 
consequence of the curvature of the Fermi surface, and 
dominates the low energy behavior of the theory at the 
critical point. The reason for this disagreement is that in 
Ref.fssjl the effects of Landau damping are ignored. We 



find that these effects are crucial. 

On the other hand, Nilsson and Castro Neto^ ap- 
proach the nematic quantum phase transition within the 
more traditional methods found in the Fermi liquid the- 
ory literature. They first find an order parameter theory 
by constructing a path integral, whose classical equations 
of motion give the coUisionless Boltzmann equation as in 
high dimensional bosonization, and integrating out the 
non-critical modes. Their results, however, agree with 
ours in all essential details, including the existence of a 
z — 3 mode. They then calculate the fermion lifetime 
using the Bethe-Salpeter equations and Fermi's Golden 
rule finding that cx e^/^, as in H3.22|l . and agree with 
us, in the perturbative regime, in concluding that this 
represents both a breakdown of Fermi liquid theory and 
perturbation theory. 
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APPENDIX A: SUMMARY OF BOSONIZATION 
IN D-DIMENSIONAL FERMI SYSTEMS 

Consider the Fermi liquid theory of spinless fermions 
interacting via a short but perhaps finite range forward 
scattering interaction living in a translationally invariant 
D-dimensional world. This is a low energy theory and 
as such, following Landau we shall linearize the energy 
dispersion near the Fermi surface, though corrections to 
this may be considered when necessary. To this end, let 
us build a construction in which we linearize within N 
equally sized patches approximating the Fermi surface. 
For this construction to be reasonable, our end result 
should be relatively insensitive to the details belonging 
to this partitioning. Keep in mind, however, that a re- 
markable property of Fermi liquids is that they only re- 
quire a few of the lowest angular momentum Fermi liquid 
parameters to understand a wide range of phenomena. 

The number of patches, N, approximating the Fermi 
surface will naturally be inversely proportional to its cur- 
vature. As such, if the density, n — > cx) then — > oo 
and the construction becomes exact. An exact solution 
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to leading order in A'' of this linearized theory is there- 
fore equivalent to the asymptotic low energy limit as dic- 
tated by the renormalization group. More specifically, 
the number of patches N and the patch width A ^ \/kpX, 
where A is the energy cutoff, must be related by the con- 
dition 2T:kp = NA^ required for the Fermi system to a 
have a finite density of states and curvature (see below). 
Consequently, the number of patches N must scale as 
N - Clearly N — > oo in the infrared limit 

A — > 0. Many of the results of this paper can be under- 
stood in light of this basic reasoning. 

Under our construction, the fcrmion annihilation op- 
erator becomes: 



c(x) = 7 ,Ck- 



k 

g»(ks+q)-x 

S qePs 



(Al) 



s 



cs(x)- 



where S labels the patch, Vs is the volume in k-space 
around the point kg and the new fermion operators, 
?/'5'(x), obey the canonical commutation relations 



{c5(x),cT,(x')} -<5s.T<5^(x-x') 



(A2) 



It is also important to note the Fourier transform nor- 
malizations within the construction: 



<5^(x - x') = 5r'(xt - x;)<5(x„ - <) 



oiqf(xt-Xt) 



D~l) 



E 



(A3) 



where we have introduced A to characterize the tangen- 
tial width of the patch (7VA is the area of the Fermi 
surface) and A as an ultraviolet cutoff about the Fermi 
surface. 

Now, in terms of our fermion operators, the linearized 
Hamiltonian is 



s 



(A4) 



where the free Hamiltonian density is 



Ho(x) - • : (Vc^5(x))c5(x) - 4(x)Vc5(x) : 

(A5) 

and the forward scattering interactions are described by 
(x)-^ J dx'i^5-T(x-x')<5%(x)<5nT(x') (A6) 

Here, the density fluctuations are defined by 

(5ns (x) = :"-5(x): = :4(x)c5(x): (A7) 



and throughout this description we have been using the 
usual normal ordering procedure for any operator O: 

■.0:\G) =0^0 = -.O-.-iGlOlG) (A8) 

where we take the filled Fermi sea as our ground state: 



iG>=n n <,s\o) 

S {qe-Ps|vsq<0} 



(A9) 



It was shown by HaldaneiS*ii, Castro Neto and 
FradkirkiSiiii^ and Houghton and MarstonAfi;^ that this 
Hamiltonian can be entirely described in terms of the 
electron density operators, (5ns(x) in the high density 
limit and it is quadratic in these operators. A fermion 
operator may then be constructed following well known 
ID bosonization techniques and so the theory can be 
solved exactly. (This represented a major step forward 
since the introduction of RPA by Bohm and Pines^.) 
Here we shall outline the proof of this solution, but 
with the traditional approach of point-splitting regular- 
ization, commonly used in the one-dimensional case (see 
for example44). 

The expectation value of the density operator n5(x) in 
the ground state of the Fermi sea is clearly divergent if 
we send the density of fermions to infinity. As a result, 
in the high density limit we are interested in, it is poorly 
defined. To control this divergence, we introduce the 
point-split operator: 



t|(x)=ct;(x + e/2)c5(x-e/2) 
= ■.h%{^):-{G\h%{^)\G) 



(AlO) 



Here, by short-distance it is meant a length scale short 
compared with the separation of all operators of interest 
but long compared with physical short length scales, i.e. 
It should be noted that physically point-splitting can be 
thought of as a means of |x| ^ \e\ ^ A~^. 

The divergent part may be computed explicitly: 

N 



(G|n^(x)|G) = -^5]5]e-'i 



qt 9n<0 



e^) r ^^e-^9.e.-H9„/A (AH) 



27r(e„ -l-iA-1) 

where we have implemented the ultraviolet cutoff. A, as 
a soft cutoff e~ I This is a highly anisotropic expres- 
sion in e. It vanishes if we first send e„ — s- 0, looking along 
a tangential direction (|ej 3> A~^ — > 0) but diverges 
if we first send et ^ looking along e„ ^ A~^ 0. 
Therefore, to capture the basic physics of the density 
operator we must choose the latter limit as the defini- 
tion of the point split operator. Keeping in mind that 
St{0) = ^F/(27r)^-i = 2nhvFN{0), where Ap is the 
{D — l)-dimensional Fermi surface area and N{0) is the 
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density of states at the Fermi surface, we obtain 

n-(x) . -!»^+Jn,(x) + ^H^(x) + . . . (A12) 

to leading order in e„ ^ A^^ and where we kept the 
expansion of : n^* (x) : to first order noticing the useful 
emergence of the free Hamiltonian density operator. 

Now that we have a controlled definition of the density 
operator, we proceed with computing its commutator: 

[ns (x), nf (x )J = Ss^rSt (xt - x^) 

X (^6 {xn - < - (e„ + e;)/2) - (5 (x„ - < + (e„ + ^ 

(A13) 

Expanding both sides of this equation in powers of e„ 
and equating like powers gives us the following result 

[(5ns(x), (5fiT(x')] = ~ihN{0)5s,T^s ■ V(5^(x - x') 

(A14) 

[(Sri5(x),H!^(x')] = -iftvs • V(57i5(x)(5^(x-x') 

(A15) 

Using these commutators, we may compute the equation 
of motion for the Heisenberg operator Sns{^, t): 

dtShsix, t)+vs ■ V(Sns(x, t)+ 

vs • V ^ / f^x'Fs_T(x - x')(571t(x', i) = (A16) 

where Fs-t — N{0)fs-T- This is the linearized coUi- 
sionless Boltzmann equation in operator form found by 
Castro Neto and Fradkin in the context of a coherent 
state formalism^?. We also notice through this deriva- 
tion that we may let 

^5(x) = TT^.Snli^) (A17) 



and obtain exactly the same answer. Hence, the Hamil- 
tonian may be expressed entirely in terms of the density 
operator Shsi'x.). 

A natural consequence of (|A14(I is that the density 
operator may be expressed in terms of a chiral boson 
field: 

5ns(x)=7V(0)v5-V^5(x) (A18) 

with a canonically conjugate momentum 

7r5(x) = -7V(0)vs-V(^s(x) (A19) 

These chiral bosons are a direct extension of the right / 
left chiral bosons in the context of ID bosonization. Fol- 
lowing this extension then, we may express the fermion 
operator as a vertex operator 

Vis(x) = r7s(xt)v/iV(0KA :e-*^-W/'^: ^ (A20) 
where rjsi^Ht) is a set of Klein factors responsible for en- 
suring that ■(/'s(x) obey the proper anti-commutation re- 
lations within the patch and on different patches. The 
relation between this fermion operator and the orig- 
inal C5(a;)-operators will be made precise in Section 
IHI Al via direct comparison of the perturbation series in 
F5_t(x — x') obtained in the bosonized theory. i/;s(x), 
therefore, is equivalent to cs(x) in the free case and in 
the interacting case, it is this operator projected onto the 
high-density subspace. Hence, this is an effective low en- 
ergy theory of a dense Fermi system and V's(x) can be 
viewed as actually creating Landau quasiparticles in the 
Fermi liquid phase. 

The canonical structure of the bosonized theory also 
allows us to directly write down a path integral formula- 
tion of the problem, including interactions described by 
a set of Landau parameters (for a derivation using co- 
herent states, see Ref.*]^). The action for the bosonized 
theory has the general form 



iV(0) 



/ (fxdt 

s •' 



NiO) 



r 



^ / d2a;d2a;'dii^5_T(x-x')vs-V(/3s(x)vT-V^T(x') 

mi 



S,T 



This action is a quadratic form in the Bose fields. Here 
we have not included higher order terms (such as those 
discussed in the body of the paper, in the context of the 
nematic instability). Such terms are generally present 
due to non-linearities in the fermion dispersion relation, 
as well as many-body effective interactions We have 
also not included "vertex operators" such as those as- 
sociated with pairing (BCS) interactionsiiS4i& We shall 



find it more convenient to work within this path-integral 
formulation when constructing a theory of the nematic 
quantum critical point. 

Before leaving our discussion of high dimensional 
bosonization, we should make a final comment on its va- 
lidity. As discussed in the main body of the paper, the 
boson theory here completely recovers the random phase 
approximation (RPA) in the long wavelength limit. This 
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should not be surprising since, in the asymptotic low en- 
ergy limit, both RPA and bosonization saturate the f-sum 
rule and are (formally) exact. This is a well known es- 
tablished property of bosonization, extensively discussed 
in the literature since the 1970's. 



as a result we may write near the nematic QCP: 



1 



2Ar(0) 



1 



1/4 



Mil - Ko{s) - 1/2. 



(B8) 



APPENDIX B: ANALYSIS OF K+ AND K," 



Here we present the partial fraction expansion of the 
Tective interactions 
quantum criticality. 



effective interactions V2 and V2 and their behavior near 



1. Partial fraction expansion of V^iqjS) 
In terms of the function Kq{s) we may write V^"^ as 



N{0) 
1 

2N{0) 



- ^x'^ + {1 - 62{q))x - 



(Bl) 



where x = Kq{s). The denominator is thus a cubic poly- 
nomial in Ko{s) and solving for the poles, x = St{q), we 
find 

. . / 1 / w ^ ^ 62{qf + l2d2{q)-12 



+ e'""/3^(,52(g))) 



€{-1,1,3} (B2) 



where 



h{x) = (I2V12 - 36a; + 42a;2 + 3a;3-a;(a;^ + 18a;+36))^/^ 

(B3) 

and we assign (3 = {a, b, c} to each of these three poles. 
The partial fraction expansion is 



(B4) 



with residues 



1 + 



^ {St - S+){6+ - d+) 

1 {1+s+r 



(B5) 



{6t - S+)iSt - S+) 

For |(52((?)| ^ 1 these formulas simplify to: 

St^t, S+^-t, St ^62/2 (B6) 

Zt^i/2, Z+^-i/2, Z+«l/2 (B7) 



2. Partial fraction expansion of (?> s) 

We may write V'2^(g, s) as 



V2-{q,s) 



7V(0) 



52{q)-Ko{s)[l-l 



l + Kois) 



{4-S2{q))N{0) 



(l + x) 



rX^-2 



S2{q) _ S2 
4 -So 4 - Jo 



(B9) 

The denominator is simply quadratic and we find poles 
^/Mq) x-t„\^ VMq) 



a,t X = 6^ with 



Saiq) 



St {<!) = - 



(BIO) 



Hence, we can now expand in partial fractions to obtain 

1 / Zail) . Zbil) 



^2 (l,s) ^^y^ _ ^^^^^ + _ ^^^^^^ 

(Bll) 

with residues 



2- = 



^-So 



a + Sa? 



4 -(5, 



Sa -S,^ J 

Again, these formulas simplify for \S2\ <^ 1: 

j-_ \/S2 s-_ y/h 

_ 1 _ 1 



a + Sb) 
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S,, - Sa 



4x/^' 



^Vs^ 



Thus, near the Nematic QCP we can write: 

1 / 1 



V2-{q,s) 



4N{0) \s2 + M?) 



(B12) 

(B13) 
(B14) 

(B15) 



APPENDIX C: INCLUSION OF Fo{q) 

Here we calculate Vs^t for the case when both Fq and 
F2 are present and we let F2 approach the Nematic QCP 
{F2 ^ —1). We shall find, however, that the critical 
behavior is utterly independent of i^o' Now, let 

Fs-Hq) = J^[po + F2{q) cos (2(^s - ^t))) (CI) 
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Returning to the cr"-theory, whose correlators are the 
RPA mteraction Vs.t, wc find in this case 



cPqdoj 



(2^)2 



a+(q,c.).(V+) ' -a+i-c^.-u) 

(C2) 



^2 k2 r 



where cr+ — (cq ,0-2 ) ^-^d 



(C3) 



As a result, is completely unaffected by the presence 
of Fo due to the fact that there is no field and opposite 
"signs" decouple. 

Again, we may write this expression as a function of 
X = Kq{s) and taking the inverse, we find 





lini{i^,u„}^o- 








0.0053 ln(!^)/;/^/^ 




— 0.0018 ln^(i/) 


-0.0018 ln^(!^) 




0.0077/i/^^^ 


—0.0077/!^^''^ 


Iu{Un) 


0.0280/m„' 


0.0140/-u,V^ 




0.0052 \n^{u„) 


-0.0052 ln^(ii„) 



TABLE L Low frequency / long wavelength limit of the re- 
quired integrals 



+ ax^ + bx + c 



W{x) 



(C4) 



where 
W{x) 



1 



2Sq + 82 



+ {52/2)x^ -{l-52)x + 52/2 x{l~x^) 



x{l-x') 



Soil + x)' 



and 



(2 - 5o)52 
2^0 + ^2-4' 



2(5o + 82 — 25q52 
2^o+'52-4 ' 



(C5) 



-^0<^2 



2(5o + (52-4 
(C6) 

with 5o = 1 + 1 /-Fq and 82 = 1 + 1 /F2 . Again we find that 
the polynomial is cubic. Hence, the addition of Fq only 
complicates the algebra, but not the general structure of 
the solution. We therefore continue as in Appendix IbI 

Solving the cubic equation in ljC4|l leads to an alge- 
braically complicated result. However, it simplifies near 
the nematic QCP and we find to lowest order in 62'. 



St 



So 



So- 2' 



So 



So -2' 



St = S2I2 (C7) 



which return to our previous result if we let -Fg ~^ or 
(5o ^ 00. Notice that, to lowest order S^ is independent 
of -Fqi as it should since the s-wave mode is non-critical 
and this "pole" precisely represents the critical behavior. 
Expanding in partial fractions leads to the expression 



V = 



1 



E 



(C8) 



with residue matrices 

/- 1 



/ I ^0-1 

i j ^ <5o (i5o-2) i5o (<5o-2 

zfer 1 + 7= 



y <5o (^i5o — 2 ) 



z+ = 



-i I So ( ^o(^So-2) 



^50(50-2) 

<5o('5o-2) 



(C9) 



1 (Q 

2 lo 1 



which again returns to the previous result if we let 
Fq 0. Notice that Z+ is independent of Fq to this 
order in 82- Consequently, the critical behavior of the 
fermions, H3.27|l . is unaffected by the inclusion of a non- 
critical mode such as Fq. 



APPENDIX D: EQUAL-POSITION BOSON 
PROPAGATOR 



Here we are interested in the quantity 
Gs(s.s)(0,<) = GO,(s^s)(0,<)+ 



(27r)3 



(Dl) 



the equal position part of Eq. H3.1|) . In particular, we fo- 
cus on the second term, labeling it lnZs(0,t). We begin 
by writing the interaction in terms of its spectral func- 
tion: 



de' B{k,e') 

2tt e — e' + ie' signe' 



(D2) 
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where 



and the residue of the double pole is 



B{k,uj) = B+{k,uj) coa^ (2(6*5 - 0)) 

+ {k, uj) sin^ (2{0s - (f>)) (D3) 

and (f) the direction of k. Note: in the appropriate scaling 
within a patch, — > 6*5 + 7r/2 so that the cosine factor 
scales to 1 while the sine factor scales to driving the 
B2(k,uj) contribution irrelevant (this has been checked 
explicitely). 

This allows us to do the e-integration immediately and 
rewrite our expression in a more physical form. The re- 
sult is (after letting e' — > 0): 

\nZs{0,\t\) = 



.A/2 


dkt 


.A/2 


dkn f 


/-A/2 


2^. 


/-A/2 


~^ Jo 



.A/2 


dkt 


r^'^ dk„ r°° 


/-A/2 


27r 


Jo 27r 



which we have written directly in terms of the patch co- 
ordinates. The residue of the single pole is: 



Ri{kn,e,\t\) = 



-ie\t\ _ ^ 



(e — vpkn + it signe) 



1 d 

vp dk„ 



e — vpkn 



le signe 



(D5) 
(D6) 



i?2(A:„,£, |t| 



1 d 

Vp dkn 



e — Vpkn + ie signe 



(D7) 



so that it is clear that an integration-by-parts removes the 
double pole all together, leaving us with an integral only 
over a single pole. One can view this as a cancelation of 
the double pole's contribution to the integral and expect 
the result to be less singular than the equal-time case. 

Using the symmetries of the spectral function, we may 
rewrite Eq. (|D4|I as 



lnZs(0,|t|) 



dkt 



dkn 



de 



B{k,s)x 



^" ^0 

i?i(/c„,e, \t\)+Ri{-kn,e, \t\)+R2{kn,e, |t|)-i?2(/c„, -e, |t| 

(D8) 



Performing the integration by parts, and a quick change 
of variables allows us to separate the time dependence in 
terms of the following integrals: 



\nZs{0,\t\) = 



2 dv 
Jo 
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r 



(D9) 



where we have sent the cutoffs to infinity, keeping the 
lowest order correction in A (having checked that higher 
orders are non-singular), and where the various integrals 
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(D13) 
(DM) 



20 



In this latest form of In Zs, we have defined u = k-y/K, 
// — Sy^/vF, t — vptj yfn, due to the form of the spectral 
function at a nematic instability: 

^I'^'^J - ^2 + ^2„2^fc6/4 - j,2+j,6/4 y^^"") 

where we have only included the effects of the over- 
damped mode since it is responsible for the leading sin- 
gular behavior. Here we have replaced the angular factor 
cos^(26'5) by a constant of order unity. (This is consis- 
tent since we will be interested in the angular- averaged 
fermion Green function and the angular factor in the ex- 
ponent vanishes only on a set of measure zero.) 

Table m shows the low frequency-long wavelength limit 
of these integrals computed numerically. Two basic forms 
emerge: a Xx^iy) divergence and a v~^l^ power-law diver- 
gence (with logarithmic corrections) each form occurring 
at different scales. However, the In [y] contributions all 



vanish. With these limits in mind, we perform the final 
integral and obtain 

Zs(0,|t|) = Cpexp|-A(l + V3*)i^p|^| (D16) 

valid for |f| 3> 1. In this expression, 

C.^exp/ ^ f0.0476+^ + ...n>l 

'^\iV(0)t;Fv^ V / J 

(D17) 

0.00724 1 
A = — oc (D18) 

7V(0)i;fkA kpXn 

= -4.95 + l.Sli (D19) 

where the numbers calculated for Gp were also computed 
numerically. 
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